Exploiting within-breed variability in the autochthonous Reggiana breed identified several candidate genes affecting pigmentation-related traits, stature and udder defects in cattle.

Autochthonous cattle breeds constitute important reservoirs of genetic diversity. Reggiana is an Italian local cattle breed reared in the north of Italy for the production of a mono-breed Parmigiano-Reggiano cheese. Reggiana cattle usually have a classical solid red coat colour and pale muzzle. As part of the strategies designed for the sustainable conservation of this genetic resource, we investigated at the genome-wise level the within-breed detected variability of three pigmentation-related traits (intensity of red coat colour, based on three classes - light/diluted, normal and dark; spotted patterns/piebaldism that sometime emerge in the breed; muzzle colour - pink/pale, grey and black), stature, presence/absence and number of supernumerary teats and teat length. A total of 1776 Reggiana cattle (about two-thirds of the extant breed population) were genotyped with the GeneSeek GGP Bovine 150k SNP array and single-marker and haplotype-based GWASs were carried out. The results indicated that two main groups of genetic factors affect the intensity of red coat colour: darkening genes (including EDN3 and a few other genes) and diluting genes (including PMEL and a few other genes). Muzzle colour was mainly determined by MC1R gene markers. Piebaldism was mainly associated with KIT gene markers. Stature was associated with BTA6 markers upstream of the NCAPG-LCORL genes. Teat defects were associated with TBX3/TBX5, MCC and LGR5 genes. Overall, the identified genomic regions not only can be directly used in selection plans in the Reggiana breed, but also contribute to clarifying the genetic mechanisms involved in determining exterior traits in cattle.

Single-marker and haplotype-based genome-wide association studies for the number of teats in two heavy pig breeds.

The number of teats is a reproductive-related trait of great economic relevance as it affects the mothering ability of the sows and thus the number of properly weaned piglets. Moreover, genetic improvement of this trait is fundamental to parallelly help the selection for increased litter size. We present the results of single-marker and haplotypes-based genome-wide association studies for the number of teats in two large cohorts of heavy pig breeds (Italian Large White and Italian Landrace) including 3990 animals genotyped with the 70K GGP Porcine BeadChip and other 1927 animals genotyped with the Illumina PorcineSNP60 BeadChip. In the Italian Large White population, genome scans identified three genome regions (SSC7, SSC10, and SSC12) that confirmed the involvement of the VRTN gene (as we previously reported) and highlighted additional loci known to affect teat counts, including the FRMD4A and HOXB1 gene regions. A different picture emerged in the Italian Landrace population, with a total of 12 genome regions in eight chromosomes (SSC3, SSC6, SSC8, SSC11, SSC13, SSC14, SSC15, and SSC16) mainly detected via the haplotype-based genome scan. The most relevant QTL was close to the ARL4C gene on SSC15. Markers in the VRTN gene region were not significant in the Italian Landrace breed. The use of both single-marker and haplotype-based genome-wide association analyses can be helpful to exploit and dissect the genome of the pigs of different populations. Overall, the obtained results supported the polygenic nature of the investigated trait and better elucidated its genetic architecture in Italian heavy pigs.

Haplotype-based genome-wide association studies reveal new loci for haematological and clinical-biochemical parameters in Large White pigs.

We report haplotype-based GWASs for 33 blood parameters measured in 843 Italian Large White pigs. In the single-trait analysis, a total of 30 QTL for number of basophils, six erythrocyte traits (haemoglobin, haematocrit, mean corpuscular haemoglobin, mean corpuscular haemoglobin concentration, mean corpuscular volume and red blood cell count) and two clinical-biochemical traits (alkaline phosphatase and Ca2+ contents) were identified. In the multiple-trait analysis, a total of five QTL affected three different clusters of traits. Only four of these QTL were already reported in the single-marker and multi-marker GWASs we previously carried out on the same pig population. QTL on SSC11 and SSC17 showed effects on multiple traits. These results further dissected the genetic architecture of parameters that could be used as proxies in breeding programmes for more complex traits. In addition, these results might help to better define the pig as an animal model for several blood-related biological functions.

Genome-wide association studies for the number of teats and teat asymmetry patterns in Large White pigs.

The number of teats is a morphological trait that influences the mothering ability of the sows and thus their reproduction performances. In this study, we carried out GWASs for the total number of teats and other 12 related parameters in 821 Italian Large White heavy pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip array. For four investigated parameters (total number of teats, the number of teats of the left line, the number of teats of the right line and the maximum number of teats comparing the two sides), significant markers were identified on SSC7, in the region of the vertnin (VRTN) gene. Significant markers for the numbers of posterior teats and the absolute difference between anterior and posterior teat numbers were consistently identified on SSC6. The most significant SNP for these parameters was an intron variant in the TOX high mobility group box family member 3 (TOX3) gene. For the other four parameters (absolute difference between the two sides; anterior teats; the ratio between the posterior and the anterior number of teats; and the absence or the presence of extra teats) only suggestively significant markers were identified on several other chromosomes. This study further supported the role of the VRTN gene region in affecting the recorded variability of the number of teats in the Italian Large White pig population and identified a genomic region potentially affecting the biological mechanisms controlling the developmental programme of morphological features in pigs.

Genome-wide association studies for iris pigmentation and heterochromia patterns in Large White pigs.

Eye colour genetics have been extensively studied in humans since the rediscovery of Mendel's laws. This trait was first interpreted using simplistic genetic models but soon it was realised that it is more complex. In this study, we analysed eye colour variability in a Large White pig population (n = 897) and report the results of GWASs based on several comparisons including pigs having four main eye colour categories (three with both pigmented eyes of different brown grades: pale, 17.9%; medium, 14.8%; and dark, 54.3%; another one with both eyes completely depigmented, 3.8%) and heterochromia patterns (heterochromia iridis - depigmented iris sectors in pigmented irises, 3.2%; heterochromia iridum - one whole eye iris of depigmented phenotype and the other eye with the iris completely pigmented, 5.9%). Pigs were genotyped with the Illumina PorcineSNP60 BeadChip and GEMMA was used for the association analyses. The results indicated that SLC45A2 (on chromosome 16, SSC16), EDNRB (SSC11) and KITLG (SSC5) affect the different grades of brown pigmentation of the eyes, the bilateral eye depigmentation defect and the heterochromia iridis defect recorded in this white pig population respectively. These genes are involved in several mechanisms affecting pigmentation. Significant associations for the eye depigmented patterns were also identified for SNPs on two SSC4 regions (including two candidate genes: NOTCH2 and PREX2) and on SSC6, SSC8 and SSC14 (including COL17A1 as candidate gene). This study provided useful information to understand eye pigmentation mechanisms, further valuing the pig as animal model to study complex phenotypes in humans.

Comparative evaluation of genomic inbreeding parameters in seven commercial and autochthonous pig breeds.

Single nucleotide polymorphism (SNP) genotyping tools, which can analyse thousands of SNPs covering the whole genome, have opened new opportunities to estimate the inbreeding level of animals directly using genome information. One of the most commonly used genomic inbreeding measures considers the proportion of the autosomal genome covered by runs of homozygosity (ROH), which are defined as continuous and uninterrupted chromosome portions showing homozygosity at all loci. In this study, we analysed the distribution of ROH in three commercial pig breeds (Italian Large White, n = 1968; Italian Duroc, n = 573; and Italian Landrace, n = 46) and four autochthonous breeds (Apulo-Calabrese, n = 90; Casertana, n = 90; Cinta Senese, n = 38; and Nero Siciliano, n = 48) raised in Italy, using SNP data generated from Illumina PorcineSNP60 BeadChip. We calculated ROH-based inbreeding coefficients (FROH) using ROH of different minimum length (1, 2, 4, 8, 16 Mbp) and compared them with several other genomic inbreeding coefficients (including the difference between observed and expected number of homozygous genotypes (FHOM)) and correlated all these genomic-based measures with the pedigree inbreeding coefficient (FPED) calculated for the pigs of some of these breeds. Autochthonous breeds had larger mean size of ROH than all three commercial breeds. FHOM was highly correlated (0.671 to 0.985) with FROH measures in all breeds. Apulo-Calabrese and Casertana had the highest FROH values considering all ROH minimum lengths (ranging from 0.273 to 0.189 and from 0.226 to 0.152, moving from ROH of minimum size of 1 Mbp (FROH1) to 16 Mbp (FROH16)), whereas the lowest FROH values were for Nero Siciliano (from 0.072 to 0.051) and Italian Large White (from 0.117 to 0.042). FROH decreased as the minimum length of ROH increased for all breeds. Italian Duroc had the highest correlations between all FROH measures and FPED (from 0.514 to 0.523) and between FHOM and FPED (0.485). Among all analysed breeds, Cinta Senese had the lowest correlation between FROH and FPED. This might be due to the imperfect measure of FPED, which, mainly in local breeds raised in extensive production systems, cannot consider a higher level of pedigree errors and a potential higher relatedness of the founder population. It appeared that ROH better captured inbreeding information in the analysed breeds and could complement pedigree-based inbreeding coefficients for the management of these genetic resources.

Comparative selection signature analyses identify genomic footprints in Reggiana cattle, the traditional breed of the Parmigiano-Reggiano cheese production system.

Reggiana is an autochthonous cattle breed reared mainly in the province of Reggio Emilia, located in the North of Italy. Reggiana cattle (originally a triple-purpose population largely diffused in the North of Italy) are characterised by a typical solid red coat colour. About 2500 cows of this breed are currently registered to its herd book. Reggiana is now considered a dual-purpose breed even if it is almost completely dedicated to the production of a mono-breed branded Protected Designation of Origin Parmigiano-Reggiano cheese, which is the main driver of the sustainable conservation of this local genetic resource. In this study, we provided the first overview of genomic footprints that characterise Reggiana and define the diversity of this local cattle breed. A total of 168 Reggiana sires (all bulls born over 35 years for which semen was available) and other 3321 sires from 3 cosmopolitan breeds (Brown, Holstein and Simmental) were genotyped with the Illumina BovineSNP50 panel. ADMIXTURE analysis suggested that Reggiana breed might have been influenced, at least in part, by the other three breeds included in this study. Selection signatures in the Reggiana genome were identified using three statistical approaches based on allele frequency differences among populations or on properties of haplotypes segregating in the populations (fixation index (FST); integrated haplotype score; cross-population extended haplotype homozygosity). We identified several regions under peculiar selection in the Reggiana breed, particularly on bovine chromosome (BTA) 6 in the KIT gene region, that is known to be involved in coat colour pattern distribution, and within the region of the LAP3, NCAPG and LCORL genes, that are associated with stature, conformation and carcass traits. Another already known region that includes the PLAG1 gene (BTA14), associated with conformation traits, showed a selection signature in the Reggiana cattle. On BTA18, a signal of selection included the MC1R gene that causes the red coat colour in cattle. Other selection sweeps were in regions, with high density of quantitative trait loci for milk production traits (on BTA20) and in several other large regions that might have contributed to shape and define the Reggiana genome (on BTA17 and BTA29). All these results, overall, indicate that the Reggiana genome might still contain several signs of its multipurpose and non-specialised utilisation, as already described for other local cattle populations, in addition to footprints derived by its ancestral origin and by its adaptation to the specialised Parmigiano-Reggiano cheese production system.

A machine learning approach for the identification of population-informative markers from high-throughput genotyping data: application to several pig breeds.

Single nucleotide polymorphisms (SNPs) able to describe population differences can be used for important applications in livestock, including breed assignment of individual animals, authentication of mono-breed products and parentage verification among several other applications. To identify the most discriminating SNPs among thousands of markers in the available commercial SNP chip tools, several methods have been used. Random forest (RF) is a machine learning technique that has been proposed for this purpose. In this study, we used RF to analyse PorcineSNP60 BeadChip array genotyping data obtained from a total of 2737 pigs of 7 Italian pig breeds (3 cosmopolitan-derived breeds: Italian Large White, Italian Duroc and Italian Landrace, and 4 autochthonous breeds: Apulo-Calabrese, Casertana, Cinta Senese and Nero Siciliano) to identify breed informative and reduced SNP panels using the mean decrease in the Gini Index and the Mean Decrease in Accuracy parameters with stability evaluation. Other reduced informative SNP panels were obtained using Delta, Fixation index and principal component analysis statistics, and their performances were compared with those obtained using the RF-defined panels using the RF classification method and its derived Out Of Bag rates and correct prediction proportions. Therefore, the performances of a total of six reduced panels were evaluated. The correct assignment of the animals to its breed was close to 100% for all tested approaches. Porcine chromosome 8 harboured the largest number of selected SNPs across all panels. Many SNPs were included in genomic regions in which previous studies identified signatures of selection or genes (e.g. ESR1, KITL and LCORL) that could contribute to explain, at least in part, phenotypically or economically relevant traits that might differentiate cosmopolitan and autochthonous pig breeds. Random forest used as preselection statistics highlighted informative SNPs that were not the same as those identified by other methods. This might be due to specific features of this machine learning methodology. It will be interesting to explore if the adaptation of RF methods for the identification of selection signature regions could be able to describe population-specific features that are not captured by other approaches.

Signatures of de-domestication in autochthonous pig breeds and of domestication in wild boar populations from MC1R and NR6A1 allele distribution.

Autochthonous pig breeds are usually reared in extensive or semi-extensive production systems that might facilitate contact with wild boars and, thus, reciprocal genetic exchanges. In this study, we analysed variants in the melanocortin 1 receptor (MC1R) gene (which cause different coat colour phenotypes) and in the nuclear receptor subfamily 6 group A member 1 (NR6A1) gene (associated with increased vertebral number) in 712 pigs of 12 local pig breeds raised in Italy (Apulo-Calabrese, Casertana, Cinta Senese, Mora Romagnola, Nero Siciliano and Sarda) and south-eastern European countries (Krskopolje from Slovenia, Black Slavonian and Turopolje from Croatia, Mangalitsa and Moravka from Serbia and East Balkan Swine from Bulgaria) and compared the data with the genetic variability at these loci investigated in 229 wild boars from populations spread in the same macro-geographic areas. None of the autochthonous pig breeds or wild boar populations were fixed for one allele at both loci. Domestic and wild-type alleles at these two genes were present in both domestic and wild populations. Findings of the distribution of MC1R alleles might be useful for tracing back the complex genetic history of autochthonous breeds. Altogether, these results indirectly demonstrate that bidirectional introgression of wild and domestic alleles is derived and affected by the human and naturally driven evolutionary forces that are shaping the Sus scrofa genome: autochthonous breeds are experiencing a sort of 'de-domestication' process, and wild resources are challenged by a 'domestication' drift. Both need to be further investigated and managed.

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool-seq approach.

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool-seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5-10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.

Genome-wide association study for the level of serum electrolytes in Italian Large White pigs.

Calcium, magnesium and phosphorus are essential electrolytes involved in a large number of biological processes. Imbalance of these minerals in blood may indicate clinically relevant conditions and are important in inferring acute or chronic pathologies in humans and animals. In this work, we carried out a genome-wide association study (GWAS) for the level of these three electrolytes in the serum of 843 performance-tested Italian Large White pigs. All pigs were genotyped with the Illumina PorcineSNP60 BeadChip, and GWAS was carried out using genome-wide efficient mixed-model association. For the level of Ca(2+) , eight single nucleotide polymorphisms (SNPs) were significant, considering a false discovery rate (FDR) < 0.05, and another eight were above the moderate association threshold (Pnominal value  < 5.00E-05). These SNPs are distributed in four porcine chromosomes (SSC): SSC8, SSC11, SSC12 and SSC13. In particular, a few putative different signals of association detected on SSC13 and one on SSC12 were in genes or close to genes involved in calcium metabolism (P2RY1, RAP2B, SLC9A9, C3orf58, TSC22D2, PLCH1 and CACNB1). Only one SNP (on SSC7) and six SNPs (on SSC2 and SSC7) showed moderate association with the level of magnesium and phosphorus respectively. The association signals for these two latter minerals might identify genes not known thus far for playing a role in their biological functions and regulations. In conclusion, our GWAS contributed to increased knowledge on the role that calcium, magnesium and phosphorus may play in the genetically determined physiological mechanisms affecting the natural variability of mineral levels in mammalian blood.

Metabolomics evidences plasma and serum biomarkers differentiating two heavy pig breeds.

In pigs, many production traits are known to vary among breeds or lines. These traits can be considered end phenotypes or external traits as they are the final results of complex biological interactions and processes whose fine biological mechanisms are still largely unknown. This study was designed to compare plasma and serum metabolomic profiles between animals of two heavy pig breeds (12 Italian Large White and 12 Italian Duroc), testing indirectly the hypothesis that different genetic backgrounds might be the determining factors of differences observed on the level of metabolites in the analyzed biofluids between breeds. We used a targeted metabolomic approach based on mass spectrometric detection of about 180 metabolites and applied a statistical validation pipeline to identify differences in the metabolomic profiles of the two heavy pig breeds. Blood samples were collected after jugulation at the slaughterhouse and prepared for metabolomics analysis that was carried out using the Biocrates AbsoluteIDQ p180 Kit, covering five different biochemical classes: glycerophospholipids, amino acids, biogenic amines, hexoses and acylcarnitines. A statistical pipeline that included the selection of the most relevant metabolites differentiating the two breeds by sparse Partial Least Squares Discriminant Analysis (sPLS-DA) was coupled with a stability test and significance test determined with leave one out and permutation procedures. sPLS-DA plots clearly separated the pigs of the two investigated breeds. A few metabolites (a total of five metabolites considering the two biofluids) involved in key metabolic pathways largely contributed to these differences between breeds. In particular, a higher level of the sphingomyelins SM (OH) C14:1 (both in plasma and serum), SM (OH) C16:1 (in serum) and SM C16:0 (in serum) were observed in Italian Duroc than in Italian Large White pigs and the inverse was for the biogenic amine kynurenine (in plasma). The level of another biogenic amine (acetylornithine) was higher in Italian Large White than in Italian Duroc pigs in both analysed biofluids. These results provided biomarkers that could be important to understand the biological differences between these two heavy pig breeds. In particular, according to the functional role played by sphingomyelins in obesity-induced inflammatory responses, it could be possible to speculate that a higher level of sphingomyelins in Italian Duroc might be related to the higher interrmuscular fat deposition of this breed compared with the Italian Large White. Additional studies will be needed to evaluate the relevance of these biomarkers for practical applications in pig breeding and nutrition.

The albinism of the feral Asinara white donkeys (Equus asinus) is determined by a missense mutation in a highly conserved position of the tyrosinase (TYR) gene deduced protein.

A feral donkey population (Equus asinus), living in the Asinara National Park (an island north-west of Sardinia, Italy), includes a unique white albino donkey subpopulation or colour morph that is a major attraction of this park. Disrupting mutations in the tyrosinase (TYR) gene are known to cause recessive albinisms in humans (oculocutaneous albinism Type 1; OCA1) and other species. In this study, we analysed the donkey TYR gene as a strong candidate to identify the causative mutation of the albinism of these donkeys. The TYR gene was sequenced from 13 donkeys (seven Asinara white albino and six coloured animals). Seven single nucleotide polymorphisms were identified. A missense mutation (c.604C>G; p.His202Asp) in a highly conserved amino acid position (even across kingdoms), which disrupts the first copper-binding site (CuA) of functional protein, was identified in the homozygous condition (G/G or D/D) in all Asinara white albino donkeys and in the albino son of a trio (the grey parents had genotype C/G or H/D), supporting the recessive mode of inheritance of this mutation. Genotyping 82 donkeys confirmed that Asinara albino donkeys had genotype G/G whereas all other coloured donkeys had genotype C/C or C/G. Across-population association between the c.604C>G genotypes and the albino coat colour was highly significant (P = 6.17E-18). The identification of the causative mutation of the albinism in the Asinara white donkeys might open new perspectives to study the dynamics of this putative deleterious allele in a feral population and to manage this interesting animal genetic resource.

Deconstructing the pig sex metabolome: Targeted metabolomics in heavy pigs revealed sexual dimorphisms in plasma biomarkers and metabolic pathways.

Metabolomics has opened new possibilities to investigate metabolic differences among animals. In this study, we applied a targeted metabolomic approach to deconstruct the pig sex metabolome as defined by castrated males and entire gilts. Plasma from 545 performance-tested Italian Large White pigs (172 castrated males and 373 females) sampled at about 160 kg live weight were analyzed for 186 metabolites using the Biocrates AbsoluteIDQ p180 Kit. After filtering, 132 metabolites (20 AA, 11 biogenic amines, 1 hexose, 13 acylcarnitines, 11 sphingomyelins, 67 phosphatidylcholines, and 9 lysophosphatidylcholines) were retained for further analyses. The multivariate approach of the sparse partial least squares discriminant analysis was applied, together with a specifically designed statistical pipeline, that included a permutation test and a 10 cross-fold validation procedure that produced stability and effect size statistics for each metabolite. Using this approach, we identified 85 biomarkers (with metabolites from all analyzed chemical families) that contributed to the differences between the 2 groups of pigs ( < 0.05 at the stability statistic test). All acylcarnitines and almost all biogenic amines were higher in castrated males than in gilts. Metabolites involved in tryptophan catabolism had the largest differences (i.e., delta = 20% for serotonin) between castrated males (higher) and gilts (lower). The level of several AA (Ala, Arg, Gly, His, Lys, Ser, Thr, and Trp) was higher in gilts (delta was from approximately 1.0 to approximately 4.8%) whereas products of AA catabolism (taurine, 2-aminoadipic acid, and methionine sulfoxide) were higher in castrated males (delta was approximately 5.0-6.0%), suggesting a metabolic shift in castrated males toward energy storage and lipid production. Similar general patterns were observed for most sphingomyelins, phosphatidylcholines, and lysophosphatidylcholines. Metabolomic pathway analysis and pathway enrichment identified several differences between the 2 sexes. This metabolomic overview opened new clues on the biochemical mechanisms underlying sexual dimorphism that, on one hand, might explain differences in terms of economic traits between castrated male pigs and entire gilts and, on the other hand, could strengthen the pig as a model to define metabolic mechanisms related to fat deposition.

Differentiation of meat from European wild boars and domestic pigs using polymorphisms in the MC1R and NR6A1 genes.

Wild boar meat cannot be easily distinguished from domestic pig meat, especially in processed products, thus it can be fraudulently substituted with cheaper domestic pork. In this study we genotyped polymorphisms in two genes (MC1R, affecting coat color and NR6A1, associated with number of vertebrae) in 293 domestic pigs of five commercial breeds, 111 wild boars sampled in Italy, and 90 in Slovenia and other Western Balkan regions. Allele and genotype frequency data were used to set up a DNA-based method to distinguish meat of wild boars and domestic pigs. Genotyping results indicated that domesticated genes were introgressed into wild boar populations. This complicated the determination of the origin of the meat and would cause a high error rate if markers of only one gene were used. The combined use of polymorphisms in the two analyzed genes substantially reduced false negative results.

A frameshift mutation in the melanophilin gene causes the dilute coat colour in rabbit (Oryctolagus cuniculus) breeds.

In rabbit, the dilute locus is determined by a recessive mutated allele (d) that causes the dilution of both eumelanic and pheomelanic pigmentations. In mice, similar phenotypes are determined by mutations in the myosin VA, Rab27a and melanophilin (MLPH) genes. In this study, we investigated the rabbit MLPH gene and showed that a mutation in this gene appears responsible for the dilute coat colour in this species. Checkered Giant F1 families segregating for black and grey (diluted or blue) coat colour were first genotyped for a complex indel in intron 1 of the MLPH gene that was completely associated with the coat colour phenotype (θ = 0.00; LOD = 4.82). Then, we sequenced 6357 bp of the MLPH gene in 18 rabbits of different coat colours, including blue animals. A total of 165 polymorphisms were identified: 137 were in non-coding regions and 28 were in coding exons. One of them was a frameshift deletion in exon 5. Genotyping the half-sib families confirmed the complete cosegregation of this mutation with the blue coat colour. The mutation was analysed in 198 rabbits of 23 breeds. All Blue Vienna and all other blue/grey/ash rabbits in other breeds (Californian, Castor Rex, Checkered Giant, English Spot, Fairy Marburg and Fairy Pearly) were homozygous for this deletion. The identification of MLPH as the responsible gene for the dilute locus in rabbit provides a natural animal model for human Griscelli syndrome type 3 and a new mutant to study the role of this gene on pigmentation.

Polymorphisms in an obesity-related gene (PCSK1) are associated with fat deposition and production traits in Italian heavy pigs.

The proprotein convertase subtilisin/kexin type 1 (PCSK1) gene encodes the prohormone convertase 1/3 enzyme that processes prohormones into functional hormones that, in turn, regulate central and peripheral energy metabolism. Mutations in the human PCSK1 gene cause severe monogenic obesity or confer risk of obesity. We herein investigated the porcine PCSK1 gene with the aim of identifying polymorphisms associated with fat deposition and production traits in Italian heavy pigs. By re-sequencing about 5.1 kb of this gene in 21 pigs of different breeds, we discovered 14 polymorphisms that were organized in nine haplotypes, clearly distributed in two clades of putative European and Asian origin. Then we re-mapped this gene on porcine chromosome 2 and analysed its expression in several tissues including gastric oxyntic mucosa of weanling pigs in which PCSK1 processes the pre-pro-ghrelin into ghrelin, which in turn is involved in the control of feed intake and energy metabolism. Association analyses between PCSK1 single-nucleotide polymorphisms (SNPs) and production, carcass and several other traits were conducted on five groups of pigs from three different experimental designs, for a total of 1221 animals. Results indicated that the analysed SNPs were associated (P < 0.01 or P < 0.05) with several traits including backfat thickness and visible intermuscular fat in Italian Duroc (ID) and growth performances in Italian Large White (ILW) and in ILW × Italian Landrace pigs. However, the effects estimated in the ILW were opposite to the effects reported in the ID pigs. Suggestive association (P < 0.10) was observed with muscle cathepsin B activity, opening, if confirmed, potential applications to reduce the excessive softness defect of the green hams that is of particular concern for the processing industry. The results obtained supported the need to further investigate the PCSK1 gene to fully exploit the value of its variability and apply this information in pig breeding programmes.

A whole genome scan for QTL affecting milk protein percentage in Italian Holstein cattle, applying selective milk DNA pooling and multiple marker mapping in a daughter design.

We report on a complete genome scan for quantitative trait loci (QTL) affecting milk protein percentage (PP) in the Italian Holstein-Friesian cattle population, applying a selective DNA pooling strategy in a daughter design. Ten Holstein-Friesian sires were chosen, and for each sire, about 200 daughters, each from the high and low tails of estimated breeding value for PP, were used to construct milk DNA pools. Sires and pools were genotyped for 181 dinucleotide microsatellites covering all cattle autosomes. Sire marker allele frequencies in the pools were obtained by shadow correction of peak height in the electropherograms. After quality control, pool data from eight sires were used for all subsequent analyses. The QTL heterozygosity estimate was lower than that of similar studies in other cattle populations. Multiple marker mapping identified 19 QTL located on 14 chromosomes (BTA1, 2, 3, 4, 5, 6, 8, 9, 12, 14, 17, 20, 23 and 27). The sires were also genotyped for seven polymorphic sites in six candidate genes (ABCG2, SPP1, casein kappa, DGAT1, GHR and PRLR) located within QTL regions of BTA6, 14 and 20 found in this study. The results confirmed or excluded the involvement of some of the analysed markers as the causative polymorphic sites of the identified QTL. The QTL identified, combined with genotype data of these candidate genes, will help to identify other quantitative trait genes and clarify the complex QTL patterns observed for a few chromosomes. Overall, the results are consistent with the Italian Holstein population having been under long-term selection for high PP.

A first comparative map of copy number variations in the sheep genome.

We carried out a cross species cattle-sheep array comparative genome hybridization experiment to identify copy number variations (CNVs) in the sheep genome analysing ewes of Italian dairy or dual-purpose breeds (Bagnolese, Comisana, Laticauda, Massese, Sarda, and Valle del Belice) using a tiling oligonucleotide array with ~385,000 probes designed on the bovine genome. We identified 135 CNV regions (CNVRs; 24 reported in more than one animal) covering ~10.5 Mb of the virtual sheep genome referred to the bovine genome (0.398%) with a mean and a median equal to 77.6 and 55.9 kb, respectively. A comparative analysis between the identified sheep CNVRs and those reported in cattle and goat genomes indicated that overlaps between sheep and both other species CNVRs are highly significant (P<0.0001), suggesting that several chromosome regions might contain recurrent interspecies CNVRs. Many sheep CNVRs include genes with important biological functions. Further studies are needed to evaluate their functional relevance.